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+^ ! ABSTRACT. An approximate method to calculate dielectric response and relaxation functions 

r^ \ for water saturated sedimentary rocks is tested for realistic threedimensional pore space images. 

The test is performed by comparing the prediction from the approximate method against the ex- 
^' ■ act solution. The approximate method is based on image analysis and local porosity theory. An 

empirical rule for the specification of the length scale in local porosity theory is advanced. The 
^ ! results from the exact solution are compared to those obtained using local porosity theory and 

OO ■ various other approximate mixing laws. The calculation based on local porosity theory is found 

00 ' 

_i. ■ to yield improved quantitative agreement with the exact result. 
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"O ■ 1. INTRODUCTION 

O 

o 



Applying a small electric field to a heterogeneous mixture of homogeneous and isotropic di- 
electrics can give rise to an effective dielectric behaviour that may differ substantially from 
that of the constituents. The effective dielectric constant of the mixture depends not only on 
the dielectric constants of the constituent materials, but also on the geometrical microstruc- 
ture of the mixture [|l], H]. It is tempting to utilize the dependence on the microstructure for 
deducing microstructural information from dielectric measurements. This motivates (at least 
partially) the continued interest in the dielectric response of waterfiUed sedimentary rocks or 
soils [^ ^ ^, ^, 0, |S|, 1^]. Such studies are important for the interpretation of geophysical or 
petrophysical borehole measurements [ [T0| , [11]]. Before microstructural information can be de- 
duced from dielectric measurements a reliable theory is needed that links the observed dielectric 
response with the desired microstructural information. Our objective in this paper is to discuss 
local porosity theory [|T^ as an approximate relation between microstructural information and 
dielectric response. We compare its prediction against classical mixing laws and against the 
exact value for the effective dielectric constant. Our discussion starts with reminding the reader 
how the microstructure enters via the microscopic equations. Next we review briefly the most 
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popular mixing laws. In Section 4 we recall the basics of local porosity theory , and advance a 
new length scale, called the percolation length. In Section 5 we present results for four different 
threedimensional sandstone samples. We give only a short summary of selected results. A more 
detailed account is in preparation [jTSj] . 



2. MICROSCOPIC EQUATIONS 

Microscopically the electric fields and potentials are governed by Maxwells equations in the 
quasistatic approximation. To be more specific let us consider from now on a two phase mix- 
ture of water and rock. The water fills the pore space of the rock. The pore space will be denoted 
as P, while the rock matrix will be denoted as M. The sets P, M G M^ are subsets of threed- 
imensional space, and their union S = P U M represents the sample. Within the quasistatic 
approximation the electrical potential t/(r) obeys the equation 

V-D(r) = 0, rG§,r^9M (1) 

D(r) = -£(r)V[/(r), rG§,r^9M (2) 

where D(r) is the electric displacement vector, 

e{r) = epXp(r) + ^MX^lr) (3) 

is the inhomogeneous dielectric function, and dWl = dF denotes the internal boundary between 
the two phases. Here 

r 1 : for r G G 

^G^""^ ~ \ : for r ^ G ^^ 

denotes the indicator function of a set G, and Sf, em are the dielectric constants (possibly fre- 
quency dependent) of the constituents (water and rock). For real rocks with typical pore sizes of 
order lOO/xm the quasistatic approximation remains valid for frequencies up to several lOOGHz. 



Differentiation of the dielectric function e{r) is not allowed at the internal boundary 
dF, and hence the equations must be supplemented with boundary conditions on the internal 
boundary. Assuming that there is no surface charge density, and denoting the unit normal to the 
surface by n we have 

limn ■ D(r + r^n) = limn ■ D(r — r^n), rGt^M (5) 

lim n X Vf/(r + r/n) = limn x Vf/(r — r/n), rG^M (6) 

where the second condition expresses continuity of the electric field component tangential to the 
surface. These equations must be supplemented with further boundary conditions representing 
the applied external potential at the sample boundaries dE. 

Given a solution f/(r; dWl) of the above equations all macroscopic properties of interest can 
in principle be calculated from it. An example would be the effective macroscopic dielectric 
constant e. It can be defined by averaging the solution t/(r; dM) of the microscopic problem 

U{r) = -6 VU{r) (7) 



where A formally represents a suitably defined ensemble or spatial average of a quantity A over 
microstructures. The corresponding measure on the space of microstructures generally depends 
on parameters ttj characterizing the statistical properties of the micro structure. In practice the 
microstructure is usually not known in sufficient detail. As a consequence it is impossible to 
calculate t/(r; dM), and one has to resort to approximate theories based on incomplete knowl- 
edge. This gives rise to the so called mixing laws. Mixing laws provide a relationship between 
the dielectric constants £i of the constituent phases i = 1, 2, ... and the effective macroscopic 
dielectric constant of the mixture e. They have the general form 

e = e{ei,e2,...;TCi,'JT2,...) (8) 

where ttj are parameters reflecting the dependence on the microstructure. These parameters 
determine which information about the microstructure can be obtained from measurements of 
the effective dielectric constant. 



3. MIXING LAWS AND RIGOROUS BOUNDS 



Let us illustrate the approximations involved in deriving mixing laws with a simple example. If 
we define a spatial average over the whole sample as 



f = J^J /(r)x,(r)rf^r (9) 

where 

\G\= jx^{r)dh (10) 

denotes the volume of a set G G M^, then we obtain by inserting eq. @ into eq. @ and 

applying eq. (^ _ 

D = -^P X.Vf/ - £M X„V[/ (11) 



If we assume x Vf/ = x Vf/ in the spirit of mean field theories, then we obtain by compari- 
son with eq. (0) the mixing law of the arithmetic average 

e = 0ep + (1 - (f))eM (12) 

where 

denotes the porosity, i.e. the volume fraction of pore space. 



Other popular mixing laws for homogeneous and isotropic systems, similar to eq. ([T2|), include 
the harmonic averages 

£= (0£-i + (l-0)e„i)"\ (14) 

the Clausius-Mossotti approximation with P as background phase 
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the Clausius-Mossotti approximation with M as background phase 

eB{(t>) = £m 1 - -— ^ =— ]=£m[ =T , (16) 

and the self-consistent effective medium approximation 

- + (1-0) ^ ^_ = (17) 



Ef + 2e ^M + 25 

which leads to a quadratic equation for e. In all of these mixing laws the porosity is all that 
is left to characterize the microstructure. Measurement of e combined with the knowledge of 
^M, ^p allows to deduce the porosity from such formulae. It is clear, however, that only the 



value for given in eq. (|T3|) is correct, and this raises the issue of how accurate the approximate 
mixing laws actually are. 

In such a situation it is useful to have rigorous upper and lower bounds for e. Such bounds can 
be derived by variational methods under various assumptions about the stochastic nature of the 
microstructure [jT^, [T3|, |T^. If the microstructure is known to be homogeneous and isotropic 



with bulk porosity 0, and if £p > Em, then 

eB(^)<e<Ec(^) (18) 

holds, where the upper and the lower bound are given by the Clausius-Mossotti formulae, eqs. 
( [161 ) and ([T5|). For Ef < em the bounds are reversed. It turns out that the simplest mixing laws, 
eq. (|12|) and (PH), violate these bounds. They will not be considered further. 

4. LOCAL POROSITY THEORY 



A fundamental drawback of the classical bounds and mixing laws is that they depend only on 
porosity as a single geometric parameter characterising the complex microstructures. Re- 
cently a new mixing law was developed which circumvents this restriction by incorporating 



fluctuations in porosity and connectivity [[T^, jTV), [18], ^^ ^ |T]]. The basic idea of the new ap 



proach, called local porosity theory, is to measure fundamental geometric observables (such as 
volume fraction, surface density, mean curvature density, Euler-characteristic or connectivity) 
within a bounded (compact) subset of the porous medium and to collect these measurements 
into various histograms. These histograms are then used in a generalization of the effective 
medium approximation to predict effective transport properties. 

Let IK(r, L) G S denote a cube of sidelength L centered at r. The set ]K(r, L) defines a mea- 
surement cell (window) inside of which local geometric properties such as porosity or specific 
internal surface are measured [[T^. The local porosity in this measurement cell ]K(r, L) is de- 



fined as 



^ |PnK(r.L)| 
^' ' |K(r,L)| ' ' 



The local porosity distribution yu(0, L) is defined as 

/i(0,L) = lV(5(0-0(r,L)) (20) 



m 

r 



fJ(M,-L + l) (21) 



where 5(x) denotes Diracs (5-distribution, and the summation runs over placements of the mea- 
surement cell. The integer 
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m = 

i=l 

is the number of placements of K(r , L) for a discretized sample (assumed to be a parallelepiped) 
with sidelengths Mj. Ideally all measurement cells should be disjoint [jT^], but in practice this 
would give very poor statistics. The support of /i(0, L) is the unit interval < < 1 for all L. 

The second geometrical ingredient for local porosity theory characterizes the connectivity of 
each measurement cell. We define 

(l : if (9 K_^ (r, L) -^ 9 K+^ (r, L) in K (r, L) n P 
Mr,L) = \ . (22) 

[ : otherwise 

where ^^ indicates that there exists a path inside the pore space of the measurement cell which 
connects the left boundary d K_^. (r, L) of K(r, L) perpendicular to the a;-direction to its right 
boundary d K+^. (r, L). Similarly we define Ay^r, L) for percolation in the y-direction, and 
A^ir, L) for the z-direction. It is possible to relate these quantities to the Euler-characteristic of 
P p2[ ] . Given eq. (|2^) we define the local percolation probability 



^/ , ^x _ Er ^x.(r, L)Ky{v, L)A,(r, L)5^^(^^^l) ^23) 

Er ^H{r,L) 

which gives the probability that a cell with local porosity percolates in all three directions. 
With these preparations the mixing law of local porosity theory reads [ fT^ ] 

^''^'^^ ~ ^ A(0, L)ii{<P, L)d(P + [ f^TV^i^ - ^(</'' L))ii{(P, L)d(P = (24) 



where Eb and €c are given in eqs. ( [T6| ) and ([T5|), /i in eq. (pO|), and A in eq. (p3l). 



The mixing law ( P4| ) is a generalization of the effective medium approximation. In fact, it 
reduces to eq. ( [T7| ) in the limit L ^ 0. In the limit L ^ 00 it also reduces to eq. ( [T7| ) 
albeit with in eq. (jTV]) replaced with A(0). In both limits the basic assumptions underlying all 
effective medium approaches become invalid. For small L the local geometries become strongly 
correlated, and this is at variance with the basic assumption of weak or no correlations. For 
large L on the other hand the assumption that the local geometry is sufficiently simple becomes 



invalid 0120. Hence we expect that formula (Q) will yield good results only for intermediate L. 



The question which L to choose has been discussed in the literature p3| , ^ pi] , P5| ]. Here we 
advance a new proposal. We suggest to use a length scale Lp, called the percolation length. It 



p{L)= / fi{<p,L)X{<f),L)d<p (25) 



is defined using the function 



which gives the total fraction of percolating cells at length L. Experience shows that p(L) fre- 
quently has a sigmoidal shape, and this has led us to define Lp as the length scale corresponding 
to the inflection point of p{L). Hence we define Lp through the condition 



(fp 



dL^ 



= (26) 



assuming that it is unique. The idea behind this definition is that at the inflection point the 
function p(L) changes most rapidly from its trivial value p(0) = at small L to its equally 
trivial value p{oo) = 1 at large L (assuming that the pore space percolates). We have observed 
that the length Lp is typically much larger than the correlation length [ pTj , p6| ]. 

We remark that there are two other important length scales associated with p{L). The first of 
these is the threshold length Lc defined by 

p{L,)=p, (27) 

where pc can be taken as the percolation threshold for the underlying lattice {pc ~ 0.248812 
for the simple cubic lattice pT]]) or as Pc = 1/3 for the effective medium approximation. This 
length scale is particularly important for network models which attempt to replace the complex 
microstructure by an effective lattice with similar statistical properties as the real sample. Of 
course L^ may not exist when (p = p{0) > p^. 

The second length scale Ls is the length at which p{L) approaches its asymptotic value p{oo) 
to a given degree of accuracy. For percolating samples p{oo) = 1. We define Ls through 

\p{oo) - p{Ls)\ < S (28) 

for small 6 > 0. The length Ls may be equated with the size of the so called "representative 
elementary volume" (REV)[^ required for representativity with respect to connectivity. Ls 
represents the scale of the averaging (smoothing) region that is needed to ensure that the fluctu- 
ating microscopic connectivity can be replaced with an averaged connectivity field defined on 
the continuum. The small parameter 5 controls the degree of smoothness. Naturally we expect 
Lp < Ls for small enough 5. 



5. RESULTS 



We have analyzed four samples of sedimentary rocks whose pore spaces were obtained by 
computer assisted microtomography. Images of two of the samples are given in Figures 1 and 3 
of [pT|]. Each data set consists of a threedimensional array of O's and I's indicating pore space 
P or matrix M. The dimensions of the array are Mi , M2 and Afs. Table I gives a synopsis of the 
characteristics of the four samples. 



TABLE I: Overview over properties of the data sets for four reservoir sandstones. 



Sample 


Description 


a 


Ml X M2 X M3 


Lp 





e 


A 


Berea 


10/xm 


128 X 128 X 128 


180/im 


0.1775 


9.827 


B 


coarse Sst20d 


30/im 


73 X 128 X 128 


420/im 


0.2470 


13.073 


C 


fine Sst6d 


lO/xm 


95 X 128 X 128 


100/im 


0.3200 


16.934 


D 


Fontainebleau 


7.5/im 


300 X 300 X 299 


225/im 


0.1355 


8.599 



Here a is the resolution, and Mj are the dimensionless sidelengths of the sample in units of a. 
The bulk porosity (p was defined in eq. dlBj ) and the length Lp in eq. 



We have solved the microscopic equations (@)-® numerically using the values ep ^ 87.74 and 
^M = 4.7 (in units of £0 = 8.854 ■ 10~^^F/m) for the dielectric constants, and calculated e from 
equation (0) using the averaging procedure defined in eq. @. The results are shown in the last 
column in Table I. 
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Figure 1: Derivative of p(L) (in /xm"^) defined in eq. (|25|). The abscissa is the size of 
measurement cells in /xm. The percolation length Lp corresponds to the position of the 

maximum. 



Next we measured fi((f), L) and A(0, L) as a function of (p and L. For the determination of A 
we used a Hoshen-Kopelman algorithm [E9p. Integrating the product according to eq. (E5|) and 



differentiating the result with respect to L we find the curves shown in Figure 1 . The locations 
of the maxima give the values of Lp tabulated in Table I. 



Finally we solve equation ( p4| ) iteratively to find the value of Elpt predicted by local poros- 
ity theory. We have plotted these values together with the predictions from the other mixing 
formulae in Figure 2. We emphasize that, contrary to spectral theories or network models, nei- 
ther the mixing laws nor the local porosity theory contains any free fitting parameters. While 
the Clausius-Mossotti predictions (upper and lower bounds) do not give good estimates the re- 
sults from the effective medium approximation and the local porosity theory are in much better 
agreement with the exact result. Note however that the effective medium values approach zero 
for infinite contrast, i.e. Sema — * for ep/^M — ^ oo. The values oisLPT on the other hand 



remain finite and are in similarly good agreement also in that limit [ |13| , [25D . 
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Figure 2: Comparison of approximate calculations for e evaluated from eq. ([T^) (upward 



triangles), eq.dlq) (downward triangles), eq.(|17[) (squares), and eq.(g4| (diamonds) with the 

exact result (circles). 



Of course the comparison with e at zero frequency is not sufficient to judge the quality of the 
approximations. Solutions of the frequency dependent complex dielectric function are in prepa- 
ration [[r3|], and the simultaneous comparison of real and imaginary parts is expected to provide 
further insight [ [T9| ] . Specifically, the solution of the frequency dependent inverse problem for 
local porosity theory is expected to yield information about porosity and connectivity fluctua- 
tions. We emphasize, however, that first the use of Lp must be better established. Currently we 
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view it as a successful empirical rule based on the available data. Further tests are necessary to 
corroborate or reject it. 
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